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ABSTRACT 



We use resistive magnetohydrodynamical simulations with the nested grid 
technique to study the formation of protoplanetary disks around protostars from 
molecular cloud cores that provide the realistic environments for planet forma- 
tion. We find that gaseous planetary-mass objects are formed in the early evolu- 
tionary phase by gravitational instability in regions that are decoupled from the 
magnetic field and surrounded by the injection points of the magnetohydrody- 
namical outflows during the formation phase of protoplanetary disks. Magnetic 
decoupling enables massive disks to form and these are subject to gravitational 
instability, even at ~ 10 AU. The frequent formation of planetary-mass objects 
in the disk suggests the possibility of constructing a hybrid planet formation sce- 
nario, where the rocky planets form later under the influence of the giant planets 
in the protoplanetary disk. 

Subject headings: magnetohydrodynamics — ISM: jets and outflows — stars: for- 
mation — stars: low-mass, brown dwarfs — stars: planetary systems: protoplan- 
etary disks — stars: planetary systems: formation 
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planets formed in the distant regions far from the central star? Almost all the dynamical 
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timescales for the various importan t processes in planet formation essentially scale with 



the Kepler rotation timescale (e.g., iKokubo &: Ida 120021 ). In the standard core-accretion 



scenario, there are severe timescale constr aints on the in situ formation of Jovian planets. 



even at 5 AU from the central star (e.g., Pollack et al. 19961: 



200d iHubickvi et al.ll2005l : iHori fc Ikoma M2008I : iLissauer et al. 



coma. Nakazawa. fc Emori 



20091 : iMachida et al.lboioh . 
Consequently, the formation of planets at greater distances (up to ~ 10 2 AU) seems to 
be impossible within an ob servationally reas onable timescale in the core-accretion scenario 
of planet formation (e.g., Ilda fc Linl I2004J ). On the other hand, the formation of giant 
planets due to gravitational instability of massive pro toplanetary disks ha s been extensively 
investigated by various authors (e.g., see review by iDurisen et al.l 120071 ). However, most 
analyses focused only on the evolution of hypothetical disks and cannot be applied directly 
to actual systems. The result of gravitational instability or absence of it does depend on how 
the protoplanetary disks are formed, but the formation of the protoplanetary disk is closely 
related to the formation of the central star. In other words, the initial conditions of planet 
formation by gravitational instability should be provided by the star formation process. 

The last decade has seen dramatic progress in our understanding of protostar formation 



(e.g., lAndre et al.l l2008h and now provides us with the opportunity to study the formation 
phase of protoplanetary disks. A highlight of the recent non-ideal magnetohydrodynamical 
(MHD) calculations of protostellar collapse from molecular cloud cores is the driving of 
outflows from the first cores and well-collimated fast jets from the protostars; these may 
be regarded as proof of the importance of various physical processes such as the ohmic 
dissipation of magnetic fields due to the low degree of ionization at higher density phases. 
In this Letter, we study the formation phase of protoplanetary disks in a self-consistent, 
non-ideal MHD system. Our protostellar collapse calculations start from a molecular cloud 
core, include (practically) all the realistic physical processes, and show how knowledge of the 
protostar formation process provides convincing evidence for the formation of giant planets 
in the early phase. 



2. Non-Ideal MHD Simulations 

The initial condition of our resistive MHD simulation is the critical Bonner- Ebert Sphere 
(an isothermal sphere at gravitational equilibrium) with a temperature of 10 K and a radius 
of 4750 AU. In order to initiate gravitational collapse, we increase the density uniformly by a 
factor of 3 to an initial central density of 3 x 10 6 cm 3 . The mass is 1.6 times the solar mass. In 
a typical simulation, we use a uniform rotation of the angular velocity of flmit — 1.1 x 10 _13 s _1 
and a uniform magnetic field strength of B m i t = 37 fi G in the initial state. The initial ratios 
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of thermal, rotational, and ma gnetic energy to the n egated gravitational energy are 0.3, 
0.005, 0.014, respectively (cf., iMachida et al.l 12008a ). We ado pt the nested-grid sc heme 
in order to increase the spatial resolution of the central region ( IMachida et al.ll2005l ); the 



numb er of nest-grids is typ ically 12 and each grid has n x x n y x n z = 128 x 128 x 32 



cells (IMachida et al.ll2006d ). As a result, spatial resolution of the innermost grid is 0.58AU 
and that of the outermost grid is 1200AU. 

To describe a realistic evolution of the magnetic field in protostar formation, we should 
take into account the non-ideal MHD effects of weakly ionized molecular gas. In general, 
the ambipolar diffusion is important in the low-density phase, but it is slow and not critical 
in the dynamically collapsing state, In the intermediate-density phase, the Hall term effect 
can p roduce a modest effect depending on the size distribution of dust grains (j War die &: Ng 
19991 ). In contrast, ohmic dissipation dominates in the high-density phase, and is shown 



to be the most efficient mech anism for the dissip ation of magnetic field in the magnetically 
supercritical cloud core (e.g.. iNakano et al.ll2002l ). Therefore, we model the non-ideal effects 



of the magnetic field by the effective resisti vity in the induction equation. We adopt the 
resistivity evolution of the fiducial model of IMachida et al.l ( 120071 ) that corresponds to the 
ionization equilibrium in standard molecular clouds. As the cloud core collapses, the degree 
of ionization decreases with increasing density. The gas becomes magnetically decoupled at 



densities of around 10 10 cm 3 but couples again when the t emperature exceeds a bout 1000 K. 



The adopted equation of state is the same as that used in IMachida et al.l (120091) that follows 
the r adiation hydrodynamical calculations of protostellar collapse (IMasunaga &: Inutsuka 
2000h . 



Results 



Figure 1 shows a typical "bird's eye-view" snapshot of our simulations. The timescale 
of the gravitational collapse (i.e., free-fall time) is a decreasing function of density; therefore, 
the dense central region shrinks faster than the less-dense surrounding regions. This property 
of gravitational collapse almost always leads to the successive decrease of mass inside the 
faster shrinking region in a run-away manner, This gravitational "run-away collapse" is 
decelerated by a gradual increase in the temperature of the central r e gion, and eventually 
a quas i- steady object called " the first core" is formed ( IL arson! Il969t IWinkler fc Newman 
1980aU bl; IMasunaga et al.lll998[ ). It consists mainly of hydrogen molecules and has a radial 
extent exceeding 10 AU. The molecular gas surrounding the first core continues to accrete 
onto it, resulting in a slow but monotonic increase in density and temperature at the center 
of the first core. If the initial angular momentum of the molecular cloud core is on the 
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Fig. 1.— Bird's eye- view of the result of non-ideal MHD simulation with nested grid 
technique, covering the evolution of the molecular cloud core to the protostar. The left panel 
shows the structure in the grid, level I = 8, where the high-density region (n = 10 10 cm -3 ; 
blue isodensity surface) and magnetic field lines are plotted. Two cocoon-like structures 
(brown) above and below the flattened core denote the zero-velocity surface inside of which 
the gas is outflowing from the center. The density contours (color and contour lines) and 
velocity vectors (thin white arrows) are projected in each wall surface. The right upper panel 
shows the structure in part of the 10th grid, where we can clearly see the central cavity in 
the outflowing region. The right lower panel (12th grid) shows the protoplanetary disk in 
the formation phase, and two newly formed planetary-mass objects in the disk. 
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order of the value suggested by observation, the resultant first core rotates significantly 
fast, and its formation corresponds to the onset of bipolar outflows driven by magnetic 



fields (lTomisakall200a IMachida et al.ll2006al Baneriee fc Pudritzll2006l : Machida et al.ll2008a 



Hennebelle fc Fromangil2008l ). In Figure 1, the outflow regions are shown by the two brown 
cocoon- like structure that correspond to the zero vertical velocity surfaces (v z = 0): gas 
inside the cocoons has a significant vertical velocity. 

Most of the angular momenta in gravitationally collapsi ng objects are remov ed by the 
Maxwell stress of the field, which is called magnetic braking (IMachida et al.l 120071 ). In addi- 
tion the outflowing gas carries away angular momentum during this phase. When the central 
temperature becomes sufficiently high (~ 2 x 10 3 K), the dissociation of hydrogen molecules 
becomes significant, providing effe ctive cooling that makes the core gravitationally un s table, 
trigger ing "the second collapse" (IMasunaga fc Inutsukal |2000| ; IMachida et al.l l2006al . 120071 . 
2008af ). 



The upper right panel of Figure 1 shows an enlarged view of the region inside the outflow- 
launching regions, where the resistivity is so significant that the magnetic fie ld is decoupled 



from t he gas. The outflow region envelopes the "dead zone" for magnetic field ( IMachida et al. 



2008aj) where magnetic braking is not operating. Therefore, the infalling gas in the first core 
maintains angular momentum and reaches the radius of the centrifugal barrier to form a 
circumstellar disk-like structure. The formation of this disk-like structure corresponds to 
the birth of the protoplanetary disk, which happens inside the outflow-launching region. 

The infalling from the envelope to the central region continues, and the radius of the 
outflow-lauching region increases over time, as does the outer boundary of the magnetic dead 
zone. Likewise, both the mass and the outer radius of the circumstellar disk increase over 
time. Eventually, the radius of the disk extends beyond the initial radius of the first core and 
the first core disappears; in other words, the first core transforms itself into a protoplanetary 
disk. 

An important consequence of the early formation of the protoplanetary disk is the 
gravitational fragmentation of the disk and the formation of planetary-mass objects, as 
shown in the lower right panel of Figure 1. The reason for the gravitational instability 
can be understood in terms of the ratio of the disk mass to the central object. Figure 2 
shows a schematic evolution of protostellar objects, in terms of mass. The vertical axis 
denotes mass (in units of solar mass) and the horizontal axis denotes time (in units of 
year). The red curve on the left-hand side depicts the mass of the fast collapsing region in 
the center of the molecular cloud core in the first collapse phase, which essentially defines 
the gravitationally unstable mass, and thus, corresponds to the Jeans mass. The decrease 
of this mass in the isothermal phase describes the run-away collapse of central region, a 
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Fig. 2.— Schematic diagram for the evolution of protostellar objects, in terms of mass. 
The vertical axis denotes mass (in units of solar mass) and the horizontal axis denotes time 
(in years). The red curve on the left-hand side depicts the mass of the fast collapsing region 
in the center of the molecular cloud core in the collapsing phase, which essentially defines 
the gravitationally unstable mass, and therefore corresponds to the Jeans mass. Note that 
the mass of the first core is much larger than the mass of the central protostar at its birth. 
The right-hand side describes the evolution in the main accretion phase, where gas in the 
envelope of the molecular cloud core accretes onto the central region and the protostar gains 
its mass. As the first core gradually changes into the protoplanetary disk, the mass of 
the protoplanetary disk remains larger than the mass of the protostar for a while. This 
configuration is gravitationally unstable and creates self-gravitating objects in the disk. The 
protostar mass increases monotonically and overwhelms the mass of the disk later in the 
accretion phase. 
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Fig. 3.— Left panel shows actual evolution of protostar mass and disk mass, where the 
horizontal axis denotes elapsed time in years after the formation of the protostar. The right 
panel shows time evolution of the gravitational instability indicator Q = nC$/ {irGY,). Q max 
is the maximum value in the disk and Q ave is the average value. The curve labeled R denotes 
the ratio of the most unstable wavelength to the disk radius. At £ ~ 800yr, R becomes 
smaller than unity so that the gravitationally unstable mode is possible in the disk, and 
planetary-mass objects are formed. 



characteristic of gravitational collapse of cooling gas. The slight increase of the central mass 
(the Jeans mass) corresponds to the formation of the first core and its gradual increase in a 
quasi-steady state. The second fall of the central mass corresponds to "the second collapse" 
that is triggered by the endothermic dissociation of hydrogen molecules. After most of the 
hydrogen molecules are dissociated, the adiabatic increase in the pressure of the atomic gas 
eventually overcomes the gravitational collapse and the second core forms. Therefore, the 
decrease of the central mass in the second collapse stage is closed by the formation of a 
protostar. We should note that the mass of the first core is much larger than the mass of 
the central protostar at its birth. 

The right-hand side of Figure 2 describes the evolution in the main accretion phase, 
where gas in the envelope of the molecular cloud core accretes onto the central region and 
the protostar gains mass. As the first core gradually changes into the protoplanetary disk, 
the mass of the protoplanetary disk remains large (and even larger than the protostar) 
for a while. This configuration is gravitationally unstable and the spiral arms are excited, 
which promotes gas accretion in the disk, but the accretion is not efficient enough to avoid 
gravitational fragmentation of the disk, creating self-gravitating objects in the disk. The 
typical mass of the formed objects is a fraction of the disk mass (i.e., smaller than 10" 2 solar 
masses) which corresponds to the range of Jovian planets to brown dwarfs. The protostar 
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mass increases monotonically and eventually overwhelms the mass of the disk. In effect 
Figure 2 describes why and how the gravitationally unstable protoplanetary disk is created. 

The left panel of Figure 3 shows the actual evolution of the protostar mass and disk mass, 
where the horizontal axis denotes elapsed time (in years after the formation of the protostar). 
The right panel shows the time evolution of the indicator of gravitational instability Q = 
^Cs/i^GT,), where k, Cs, G, and £ are epicyclic frequency, sound speed, gravitational 
constant, and surface density, respectively. Qmax is the maximum value in the disk and Q me 
is the average value. During the first several hundred years, Q remains smaller than unity 
and the disk satisfies the instability criterion obtained by local linear analysis, although 
the size of the disk remains too small for the unstable mode to appear, as shown by the 
curve labeled R that denotes the ratio of the most unstable wavelength to the disk radius. 
At t ~ 800 yr, R becomes smaller than unity so that the gravitationally unstable mode 
becomes feasible in the disk; this epoch does indeed correspond to the fragmentation of the 
disk. Note that the fragmentation has occurred in the region where the effective equation 
of state is almost adiabatic. If additional radiative cooling operates, it will become even 
more unstable, and thus, the fragmentation is definitely expected. In other words, the 
gravitational fragmentation shown in this early formation phase of the massive disk does not 
require efficient radiative cooling that is supposed to be importa nt in the later phase of the 



evolution w here the disk is less-massive and less- unstable (e.g., iGammie 1 12002 



Rice et al. 



20031 . 120051 ). Here we emphasize that we should expect the growth of gravitational instability 
in the disks that excessively satisfy the local instability criterion (Q < 1) for the mode 
whose wavelength is sufficiently smaller than the disk. Nevertheless, the consequence of the 
growth and the resultant fragmentation should be studied in more detail including the effect 
of irradiation from the central star that is not taken into account in the present calculations. 

The formation of planetary-mass objects and the excited spiral arm structure provide 
an efficient transfer of angular momentum in the disk and promote mass accretion onto 
the protostar, as shown by the rapid increase in the protostar mass in the left panel. This 
property will be further investigated in our next paper. 



4. Discussion 

Cooling Efficiency Since we are not solving realistic radiative transfer equation in our 
dynamical simulations, one might wonder whether each fragment can sufficiently shrink by 
radiating an excess energy obtained by compressional heating or not. Here we estimate the 
amount of energy that should be radiated away from a fragment and compare with the cooling 
capability of the fragments. In order for the self-gravity to dominate over gas pressure the 
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effective ratio of specific heats (7 e ff) should be less than 4/3, while th e ratio of specific heats 



(7) fo r the molecular gas in adiabatic evolution is mostly 7/5 (e.g., iMasunaga &: Inutsuka 



20001 ) . Thus, we can roughly calculate the amount of energy AE that should be radiated 



during the compressional motion of the planet with mass M p from the density pa, pressure 
Pd, and temperature Td in the gas disk to the average density p p and pressure P p of the 
planet as 

.Pd) \Pd 



AE 



M p fP d ^ r/ ' ^' ' " ' 



(1) 



7-1 VPp. 

If we take the density enhancement factor of 10 5 (e.g., pd = 10~ n g cm -3 and p p = 10~ 6 g cm -3 ) 
the second term is negligible compared to the first term on the right-hand side of the above 
equation. Thus, the expected time-averaged luminosity of the fragment that cools with the 
timescale At can be calculated as, 

where we assumed that the mean molecular weight is 2.3. On the other hand, the luminosity 
L p of the planet with radius R p and the surface temperature T p is 

L P = 4^ SB T< = 1.9 x 10- (jyV (3) 

where cxsb is Stephan-Boltzmann constant. Since (T/T d ) = (p/p d ) 7cff_1 , compression by a 
factor of p p /pd = 10 5 with 7 cff = 4/3 results in the increase of temperature by a factor of 
about 46. Thus, the fragment can radiate efficiently enough to collapse gravitationally by a 
factor of at least 10 5 in density within 100 years (i.e., the resultant surface temperature of the 
planet will be somewhat smaller than 10 3 K). Thus, the objects formed by the gravitational 
instability of the disk as in our simulation are expected to be dense enough to avoid being 
washed out in the dynamical environment. 



Further Evolution What is the fate of these gravitationally formed planetary-mass ob- 
jects? Does gravitational interaction with the gaseous disk lead to the migration of these 
objects? In fact, some of our simulations showed the objects falling into the central star, 
while in other cases they migrated outward but remained in the disk. Their chaotic behavior 
is not surprising, since eccentricities of the orbits of planetary-mass objects in the disk are, 
in general, not small. Their large eccentricities are due to the fact that they absorb infalling 
gas. 

In addition, the gravitational fragmentation and formation of planetary-mass objects 
repeat many times during the main accretion phase, where the gas accretion from the en- 
velope of the molecular cloud core continues. Some of the planetary-mass objects fall onto 
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the central star while some remain in the disk. Therefore, a statistical method seems to be 
required to predict the distribution of objects that survive the gas dispersal, which will be 
the scope of our subsequent paper. 



Implication for Observations Astronomical observation of the early formation phase 
of protoplanetary disks would enable us to directly test these expectation. Although our 
numerical simulations enable us to frequently observe the gravitational fragmentation of 
protoplanetary disks, in reality this process is hidden in dense gas and may remain invisible 
in optical and near-infrared wavelengths. However, radio interferometers are potentially 
capable of observing the interiors of protostellar cloud cores. The Atacama Large Millimeter 
Array (ALMA) would be an ideal tool for this type of observation. In order to directly 
observe the process described in this Letter, it is necessary to observe the early evolutionary 
phase within about 1000 years from the birth of a protostar. As the average lifetime of a 
protostar is estimated to be 10 5 year, the required time span for this observation is only 
1 % of the protostar's lifetime, meaning that we should observe, on average, about 100 
protostars to find one protostar at an early evolutionary phase. In our simulations, we almost 
always observe very time-dependent mass accretion through the gravitationally unstable disk. 
Therefore, another possible observational signature is the time variability of the protostellar 
source. We hope that the gravitational fragmentation reported in this Letter will be analyzed 
observationally and a realistic planet formation scenario will be refined in the near future. 

The frequent formation of planetary- mass objects in the disks during the formation 
phase suggests that we should consider their influence on the evolution of the protoplanetary 
disk. In particular, it suggests a new possibility for the core-accretion scenario in the region 
outside the orbit of the planetary-mass object where the shepherdi ng effect from a massiv e 



object prevents dust grains from raining out onto the central star (IMuto fc Inutsukall2009l ). 
which corresponds to the far right-side of Figure 2. Therefore, it might be interesting to 
search for a hybrid scenario of planet formation where the rocky planets form after the 
gravitational formation of the gaseous giant planets. 

This work was supported in part by Grants-in-Aid for Scientific Research from the 
MEXT of Japan (15740118, 16077202, 18540238, 18740104, 20540238, and 21740136). 
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